SE 012 698 



DOCUMENT RESUME 



ED 058 034 



AUTHOR 

TITEE 



INSTITUTI ON 
SPONS AGENCY 
PUB DATE 
NOTE 



Merrill, John R, 

[use of Computers in Introductory Physics 
Teaching. ] 

Dartmouth Coll-, Hanover, N, H- Dept, of Physics. 
National Science Foundation, Washington, D- C- 
[ 71 ] 

80p. 



EDRS PRICE 
DESCRIPTORS 



IDENTIFIERS 



MF-$0-65 HC-J3.29 

College Science; *Computer Assisted instruction; 
^Computer Oriented Programs; Educational Technology; 
♦Instructional Materials; *Physics; ♦Teaching 
Procedures 
Project COEXIST 



ABSTRACT 



This paper presents some of the preliminary results 
of Project COEXIST at Dartmouth College, an NSF sponsored project to 
investigate ways to use computers In introductory physics and 
mathematics teaching. Students use the computer in a number of ways 
on homework, on individual projects, and in the laboratory. Students 
write their own programs, whenever possible, on occasion, it is more 
useful to have the student programs written and saved by the teacher. 
Five areas of work are described; (1) illustrations from a broad 
range of study areas are introduced to show the versatility of 
computer-connected graphic display devices such as X-Y plotters and 
cathode ray terminals; (2) a simple, intuitive approach to the 
propagation of waves in dispersive media is shown; <3) a program is 
described which can be used to calculate a large number of 
interference and diffractive patterns; <4) general methods are given 
to find a map field patterns for general distribution of charges and 
current loops; and (5) solutions are found to some simple Schrodinger 
equations using a computer program- {Author/'TS) 
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ABSTRACT 

This note presents a simple, intuitive approach to the propagation of waves 
in dispersive media. The method Fourier analyzes a wave pulse, propagates each 
component according to its wave-vector dependent phase velocity, and then resyn- 
thesizes the wave train, Tlie method uses a computer to perform the algebra; a 
program which minimizes central processor time while maximizing the intuitive 
content of the calculation is given in the Appendix, 





- 2 « 



I. Introduction 

Wave propagation is usually limited to iion-dispersive media in introductory 
physics. This limitation is primarily due to the complexity of the algebra of 
dispersive media* Similarly^ Fourier analysis is usually just mentioned in 
most introductory physics courses. At most, Fourier analysis is applied to a 
few simple periodic functions. Again the limitation is essentially due to the 
complexity of the algebra involved not the sophistication of the concepts. 

Both dispersive wave theory and Fourier analysis are fairly simple in concept 
but algebraically complicated in practice. 

This note demonstrates a simple, straight-forward way to introduce wave 
propagations through dispersive media by means of Fourier transforms. A simple 
computer program is introduced which Fourier analyzes the pulse, propagates each 
wave-vector component separately and then resynthesizes the wave. The student 
supplies the initial pulse shape and a fozmiula for the wave-vector dependent 
phase velocity of the wave. The program has been used in second semester physics 
courses with students ranging from pre-medical majors through physics majors. 
Sophisticated students are able to understand the program and are able to write 
their own versions. 

Use of this computer approach has deepened the student understanding of 
non- dispersive wave theory as an approximation to real physical situations. The 
student’s intuitive understanding of expansion processes such as Fourier analysis 
has also been strengthened. 

II. Results 

Figure I shows a triangular pulse propagated through a medium of lightly 
dispersive properties. The phase velocity of the wave is slightly wave-vector 
dependant in the form VCk) ^ Vq (l+.Olk), The top curve in Figure I (curve a) 
is a plot of the initial pulse shape at t=0. The origin of the spatial 
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coordinate (X=0) lies at the center of the figure and the right and left sides 
of the figure correspond to +4IT. All the figures are scaled this way* and, 
within each figure ^ curves are offset vertically for clarity. Curve b of 
Figure I is the real part of the Fourier spectrum of the initial pulse; for the 
simple pulses shown as illustrations, there are no imaginary parts to the Fourier 
transforns. Curves c, d, e, f, and g in Figure I show the pulse after 2, 4, 6, 8 
and 10 units of time has passed. The pulse moves to the right as expected, and, 
due to the disperi .on of the propagation medium, the pulse distorts as time pro- 
gresses, Components of the wave with wave-vectors greater than the average k 
vector fk-0 in this case) propagate faster than the average of the pulse; lower 
wave-vector components fall behind. The propagation accentuates the average 
wave-vector component and produces ripples preceding and following the pulse. 
This simple model of a lightly dispersive medium corresponds roughly to the 
edges of an optical absorptioTi band. Hie present calculation does not include 
absorption itself, but the program is easily modified to include such effects. 

Figure 2 shows the same triangular pulse undergoing ten times stronger 
dispersion. Now V(k) = V^Cl+-lk) and the effects of dispersion are quite 
striking. Curves a and b are again the initial pulse and the real part of its 
Fourier transform. Curves c - f show the pulse after 2, 4, 6 and 8 units of 
time have passed. The distortion of the pulse follows the lines discussod in 
Figure I. Again the k-0 component is accentuated and higher and lower k com- 
ponents run ahead and fall behind respectively. 

Figure 5 shoWs the effects of the stronger dispersion on a simple model 
of a wave packet. Hie original pulse (curve a) is two an^ one-half cycles of a 
cosine wave. The (real p'art of the). Fourier transform^ demonstrates the peaking 
hear wave-vector con^onents of the basic cosine wave. Hie Fourier transform also 
illustrates which wave vector components must be included in order to product a 
simple pulse from a cosine wave train. Curves c, d and e (part B of Figure 3) 
show the real (or cosine) part of the pulse after propagation times of 2, 4 and 
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6 time units. Again the student observes the spreading of the pulse under the 
dispersion V(k) = V^Cl+-lk). The spreading now accentuates the wave^vector com^ 
ponents in the pulse near those of the basic cosine wave train. Figure 3C C^urves 
£ through h) shows the imaginary (or sine) part of the resynthesized pulse after 
2^ 4 and 6 time units. Figure 3 illustrates that dispersion not only effects the 
pulse envelope but also the phase relations within the pulse. The student can a 
also examine /Re^+ Im^ or the relative phase as function of position for various 
times. He can demonstrate energy conservation by examining the total energy in 
the pulse as a function of propagation time* 

Figure 4 shows the effect of propagating phonon wave packets with several 
basic wave-vectors through a one-dimensional lattice- The dispersion relation 
0 ) vs. k for lattice waves determining the wave-^vector dependent phase velocity: 
VCk)ecsin (ka/2)/k. Figure 4A shows the initial pulse (curve a) and the Fourier 
transform (curve b) for a wave packet whose basic wavelength is considerably 
longer than the lattice spacing (k^a ^ 1 so that 2na) . The Fourier trans- 
form is again peaked near k - all but a small fraction of the wave- 

vector components in the pulse have essentially the same phase velocity. The 
sharp changes in the initial pulse (i.e. the high wave-vector components) are 
influenced by dispersion, but the pulse shape remains approximately the same as 
time progresses. Curves c and d of Figure 4A show the pulse after. 12 and 24 
time units of propagation. 

Figure 4B shows the same information as Figure 4A except that the basic 
wave length of the pulse is now somewhat closer to the lattice spacing 
(k^a = 2 Now a non-negligibla fraction of the Fourier components are prop- 
agated with different velocities. The propagated pulse (curves g and h) is 
Subs t salt rally distorted, and the motion of the center of the pulse (the group 
velocity) is noticably less. 
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Figure 4C shows what happens when the basic wave-length is very near the 
lattice spacing Ck^a=3) . The Fourier transform (curve j) now has peaks near 
the Bragg reflection planes, and the propagated pulses (curves k, 1 and m having 
propagation times of 12, 24, 36 units) are highly distorted » The group velocity 
is close to zero (as suggested by dw/dk at this point) * The accentuation of 
the basic wave'-vector is apparent. 

Similar dispersive effects in other areas of ph)sics can also be examined. 
The student can observe the distortions of quantum mechanical wave packets as 
time progresses or of electromagnetic waves in plasmas and the inonosphere. An 
interesting example is the propagation of a helicon or whistler wave pulse. 

The dispersion relation for that magneto-plasma wave has proportional to k^ , 

The d. c. component of the wave packet sits still and the two dominant k values 
(±k^) in the pulse propagate away. 

Figure 5 is a simple test of* the procedure which occurs to almost all 

students. The dispersion of the medium for Figure S is zero, that is V(k) = 

The initial pulse keeps its shape when propagated eithex" forward in time 

(curves c and d for +2 and +4 time units respectively) or backward in time 

(curve e for -6 time units) , The imaginary or sine part of the resynthesized 

pulse is everywhere zero on the scale of this figure. Because the method of 

numerical integration is approximate, ’’ghost" (illegitimate replica) of the 

pulse occur for sufficiently large positive and negative wave-vectors. The 

infinite integral for the Fourier transform is approximated in the present 

program by a finite sum from -511 to +sn in steps of n/16. This approximation 

is sufficient for most , applications ; a better approximation costs more central 

1 

processor time to accomplish . 
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A few comments on the computer method may be useful- The computer strategy 
is very simple; a form of the program in the Dartmouth form of the language BASIC 
is given in the Appendix. Tlie program calculates the Fourier integral by the 
trapezoidal rule since the trapezoidal rule is more accurate than Simpson’s rule 
for nearly periodic functions"^. The limits on the integral sum and the number 
of partitions are chosen to follow reasonable detail in the pulse structure and to 
keep ’’ghosts” away from the region of intersect- To use the program, the stu- 
dent defines the original pulse shape as two functions FNR and FNI , the real and 
imaginary parts of the initial pulse- These functions are used to fill two 
vectors 3 and I( )) whose elements are called repeatedly in the Fourier 

transform routine. The use of vectors is considerably faster in central processor 
ti ...0 than repeated calls of the defined functions themselves- In programs which 
calculate Fourier integrals, every shortening of time is advantageous. 

Fourier transforms need sines and cosines repeatedly. SIN and COS functions 
are quite slow on most machines and so vectors CSC ) and CC 3) are used. The 
sines and cosines obtained this way are only accurate to one degree, but this 
accuracy is usually sufficient. The INT function used to find the right sine 
or cosine element is usually a very fast function. Another fast foCTi of sine 
and cosine, this one using the angle sum formulas, is used to fill the sine 
and cosine vectors initially. 

After calculating and storing the real and imaginary parts of the Fourier 
transform, the transform is printed out as tables. Each Fourier component is 
then propagated the correct amount using the phase velocity VCk) and the user 
supplied time. Finally, the resultant, resynthesized pulse is printed. TTie 
final result is the pulse as propagated the given time through a medium of 
supplied dispersion. On the Dartmouth GE635 system, each set of Fourier trans- 
forms C'^hi eh produces a complete spectrum or curve) takes about five central 
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processor seconds so the student can examine several different propagation times 
during one RUN . 

■Diis computer method is fairly fast and introduces the student to program- 
ming methods which conserve computer time. The student is also introduced to 
real problems in numerical integration, to the practice of Fourier transform 
theory and to the effects of wave propagation in dispersive media. 

IV, Conclusion 

This note has presented a simple computer approach to the theory of waves 
propagated through dispersive media. Hie introductory student sees dispersive 
media as an extension of and a comparison to the non-dispersive wave propagation 
usually introduced in text books. The program given in the Appendix has been 
used by numbers of students ranging from pre-meds to physics majors at the in- 
troductory . level . The response has been favorable. The program and illustrations 
of dispersive wave phenomena are useful from introductoiy physics courses 
through upper level solid state, plasma, and quantum mechanics courses. 
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FO OTNOTES 

Using the trapezoidal rule of integration, either larger limits on the integral 
sum or a finer mesh leads to more central processox time. Other integration 
schemes can be used. Gauss ian-Hermite quadrature using order one^hundred gives 
about as good results as; the trapezoidal rule] clever use of the Gaussian scheme 
seems to gain about a factor of two in CPU time. The author has not used 
Gaussian quadrature with classes since the integrations in the program are then 
complicated to describe- For discussions of more sophisticated integration 
schemes see A, H. Stroud and D, Secrest, Gaussim Quadrature Formulas , Prentice- 
H^ll, Englewood Cliffs 3 New Jersey, 1966- 

2 

For comments on this and other points having to do with numerical analysis 
see B- Carnahan, H-A- Luther, and J- 0- Wilkes, Applied Numer ical Methods, 

Wiley and Sons, Inc-, New York, 1969, 
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CAPTIONS 

Figure I. Propagation of a triangular pulse through a lightly dispersive medium, 

a) shows the original pulse j b) shows the real part of the Fourier transform of 
the initial pulse; c) through g) show the pulse after 2, 4, 6, 8, and 10 units 
of time have passed. The distortion of the pulse is due to the dispersive nature 
of the propagation mediiim. The zero of spatial and Fourier transform coordinates 
is at the center of the figure; the left and right sides of the figure correspond 
to -4II and +411 respectively. All the figures are scaled this way. Curves are 
separated vertically for clarity. 

Figure 2. Propagation of a triangular pulse through a strongly dispersive medium, 

a) shows the original pulse; b) shows the real part of the Fourier transform; 
c) through f) show the pulse after propagating 2, 4, 6, and 8 units of time. 

Fisure 3. Propagation of a wave packet through a strongly dispersive medium. 

A) shows the original pulse and the real part of its Fourier transform. The 
Fourier spectrum is peaked near plus and minus the wave-vector of the cosine 
part of the pulse shape. B) Tlis real or cosine part of the pulse after prop- 
agating 2, 4, and 6 units of time. C) The imaginary or sine part of the pulse 
after propagating 2, 4, and 6 units of time. 

Figure 4. Propagation of a phonon wave packet through a one-dimensional lattice. 

A) Shows the initial pulse, its Fourier spectrum, and its propagation for 12 
^d 24 time units. This pulse has a wavelength substantially longer. than the 
lattice spacing so the distortion is minimal. B) Shows the effects on a pulse 

whose basic wavelength is nearer the lattice spacing. The distortion is more 

noticeable and the group velocity is smaller. C) Shows the effects on a pulse 

whose basic wavelength is nearly the lattice spacing. The distortion is large, 

and the group velocity is nearly zero. 

Figure 5. Propagation of a triangular pulse through a non-dispersive medium. 

The origina.1 pulse and its Fourier transform are shown. Hie pulse after 2, 4 
and -6 units of time are also shown. The similarity in pulse shapes gives the 
student confidence in the simulation method. 
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FIGURE 1 
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FIGURE 2 
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FIGURE 3A 






FIGURE 3B 
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FIGURE 3C 
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FIGURE 4 A 
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FIGURE 4B 
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FIGURE 4C 
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FIGURE S 
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APPENDIX 



303 PRINT 




21 0 

320 REM FiNRj FNIsRE, IM PARTS OF FN. TO BE TRANSFORMED 
330 DEF FNRCX) 

340 LET FNas0 
350 LET X2sl 

360 IF CX-X2)*CX+X2>>0 THEN 280 
270 LET FNRsl-A3SCXFX2> 

280 FNEND 

390 DEF FNICX)=0 

333 

310 REM DEFINE SEVERAL PARAMETERS RELATED TO PI 

330 LET P9s3 •1415936S 

330 LET P8=3*P9 

340 LET P7SP9/16 

350 LET P6aSQRCP8> 

360 LET P5^3*P9 
370 

380 REM FILL FUNCTION VECTORS 
390 DIM RC500)«I<S00> 

400 LET 10=0 

410 FOR X0»-PS TO PS STEP P7 

420 LET 13^10+1 

430 LET R< I0)=FNR<X0> 

440 LET 1C I0>=FNICK0> 

450 NEXT X0 
460 

470 REM FILL SIN & G'OS LISTS 
430 DIM CC360>«SC360> 

490 LET P0=P8f 360 
500 LET GlsGOSCP0> 

510 LET S1=S1N€P0> 

520 LET GC0)=1 

530 LET SC0Ja0 

540 LET C2»l 

550 FOR 1=1 TO 360 

560 LET G3=G2 

570 LET C2=C2*C1 -S2n«S1 

580 LET S3=Sa*Cl-»-C3^Sl 

590 LET CCI>»C2 
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VAVES (CONTIMUEDJ 



60^ LET SfT>-Sa 
61 NEXT I 

630 

630 REM CALC. FOUR. TRANSF. FOR RANGE OF K* S 
6A0 PRINr '*K"j*'RE<F.T« )">**IM<F.T. )” 

650 DIM F<500).GC500) 

660 LET 19=0 'COUNTER 

673 FOR K0=-P5 TO P5 STEP P7 

680 LET 19=19+1 

690 LET 11=0 'REAL F.T. 

700 LET 1230 »IMAG F.T. 

713 LET 10=0 

720 FOR X03-PS TO P5 STEP P? 

730 LET 10=10+1 
740 LET K9=K0*X0 

750 LET K83lNTCCK9-INT€K9/'P@)*P0)/P0) 

760 LET C0=CCK8) 

773 LET S0=SCK8> 

783 LET R8=RC10) 

793 LET R9=ICI0) 

800 LET R1 =G0+R8-S0*R9 
810 LET Ra=C0*R9+S0*R8 
820 LET li=Il+Rl 
830 LET I2 = ia+Ra 

840 IF CX0-P5>+<X0+P5><>0 THEN 870 

850 LET Ilsll-Rl/2 

860 LET I2 = ia-R9/2 

870 NEXT X0 

880 LET I1=I1*P7/P6 

890 LET ia = I2*P7/’P6 

900 LET F(I9> = n 

910 LET G<I9)»ia 

915 IF INTC C 19*1 19-'|l >/S THEN 930 

920 PRINT K0*II,ia 
900 NEXT K0 
940 

950 REM RESYNTHESIS AFTER PROPAGATING EACH FOUR. 
960 REM ay AN AMOUNT VCKJ+T 
970 PRINT 

980 PRINT " TIME"1 
990 INPUT T1 

995 PRINT ••X"»"WAVECC0S)’»#"tfAVEISIN5" 

1000 FOR Xl=-5 TO S STEP *5 
10 T0 LET 11=0 
1080 LET ia=0 
1030 LET I0S0 

1040 FOR K0 3 -PS TO PS STEP P7 
1050 LET 10=10+1 

1060V REM >ii+*PUT VC K) BETWEEN HERE 
1370 LET V0«1+.01«K0 





COMPONENT 



- 21 - 



O 

ERIC 



WAVES 



-3- 



(CONTINUED) 



1060 
1090 
1100 
1110 
1130 
1130 
1140 
1150 
1160 
11 70 
1180 
11 90 
1200 
1210 
1820 
1830 
1840 
1853 
1260 
1870 
1260 



LET 
LET 
LET 
LET 
LET 
LET 
LET 
LET 
LET 
LET 
IF 



***AND HERE 
K9=K0*Ori -V0*T1 ) 

cl=C C K8 ; » *P8 > /P0 ) 

S0=SCK8> 

RS^Fcio ) 
a9=Q(I0) 

R 1 -C 0 *R8 +S 0 «R9 

Ra=C0*R9-S0*R8 

Il^Il+Rl 

I8sI8+Rg 

CK0-P5)*CK0+PS)<>0 THEN .1220 



-Rl/2 

•Ra/a 



LET II all 

let la^ia 

NE?CT K0 
LET HaH*P7^pQ 
LET I8aI8*P7/'P6 
PRINT Xl*Il,Ig 
NEXT XI 
OOTO 983 
END 
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ABSTRACT 

This note presents a simple program which can be used to calculate a 
large number of interference and diffraction patterns. Several applications 
of the method are discussed; radio antennas are emphasized. In class use, 
the student writes his own program and calculates intensity patterns for 
interference systems of interest to him. 
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Introduction 

This note describes some uses of the computer in practical applications 
of wave interference. Although the applications discussed are complicated 
to handle in closed, analytic form, they are easy to understand with the 
computer. The computer program propagates waves out from any number of 
wave sources to some point on a large sphere. The program then adds the 
amplitudes of all the waves (including phase) at that point, squares the 
result to get the total intensity, and averages over a complete cycle of the 
sources. One straightforward practical application of this method is to 
"broadside arrays" of radio antennas.^ 

Most students have seen broadside arrays. A radio station will place 
a mmber of antennas in an equally spaced row, the plane of the antennas 
facing a city. The station feeds the antennas with radio frequency currents 
of certain relative amplitudes and phases. The result is a highly directional 
transmission of the station's power, usually directly toward the city. 

The simplest broadside array feeds equal amplitude, in-phase currents 
the N antennnas , The antennas are often one-half wavelength apart. The 
transmission pattern is the N slit diffraction pattern. It is common to 
express the transmission pattern as a polar plot of intensity versus angle, 
around the sources. Plotting the results in this way affords the student 
a comparison to typical diffraction patterns in physics which show the 
iJ^Tensity pattern across a screen placed parallel to>.the plane of wave 
sources (or slits). 



Results 



Figure 1 shows the polar plot of the transmission pattern for two 
antennas separated by two wavelengths and fed with equal amplitude, in-phase 
currents. There are two sets of three maximum Intensity lobes in the for- 
ward and backward directions j there are broad maxima in the side directions. 
This is also a polar plot for the intensity pattern of a double slit dif- 
fraction experiment (narrow slits placed X/ 2 apart). 
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Figure 2A shows a more common broadside array = four antennas placed 
half wavelengths apart and fed equal, in-phase currents* The increased 
directionality is clear from the polar plot of intensity. Figure 2B is 
the polar plot for a **binomial broadside array” in which the currents fed 
to the antennas are in the r4.^tio of binomial coefficients.^ The trans- 
mission pattern has no power wasted in side lobes* all the intensity is 
radiated in the forward and backward directions. 

More complicated arrays are sometimes used. The "optimal array” has 
a complicated set of relative amplitudes of the feed currents « Two dimen- 
sional arrays are sometimes used, usually at high radio frequencies. The 
student can even consider Yagi-Uda arrays such as those used in TV recep- 
tion antennas. The reception pattern of an antenna system is simply re- 
lated to its transmission pattern. The method presented here can ^Iso ^ 
be applied to interferrometric radio telescopes, interference patterns 
through thin films, or the diffraction pattern of the sound emitted from a 
loud speaker. Using secondary waves, the student can model intensity 
patterns in auditoria. In all these problems, the method adds up waves 
from a finite (but sometimes large) number of distinct sources having 
different relative amplitudes and phases. 

Program 

All the cases mentioned above can be treated by a very simple com- 
puter program. Figure 3 is such a program in BASIC. The relative ampli- 
tudes of the N wave sources are stored in the vector A( ) ; the relative 
phases are strored in ); and the positions of the sources are stored 
in X( ) and Y( )• The wavelength of the waves is L, The program cal- 
culates the intensity 10 (averaged over a period of the wave) at 10 degree 
intervals (angle T) around a circle of radius 10. The average over one 
cycle in phase is performed by the FOR-NEXT loop on PO. The fundamental 
loop. which adds the wave Miplitudes from the N sources is lines 230 through 
290 (7 statements). As written, the program calculates the intensity as 
a function of angle for two source interference. Figures 1 and 2 were 
plotted directly from output of a program such as that shown in Figure 3. 
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Conclusion 

The program discussed in this note can be used with the smallest 
computers. The method is easy to explain even to the most mathematically 
unsophisticated students. The method of the program is entirely general 
and covers a very large number of applications of interference and dif- 
fraction phenomena. 
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Captions 
Figure 1 « 



Figure 2. 



Figure 3, 







Polar plot of the intensity pattern for two sources, two 
wavblengths apart. The program allows any number of sources 
to tie placed anywhere in the plane; the program then cal- 
culates d polar intensity pattern^ 

Polar plot of the intensity pattern for four element radio 
antenna broadside arrays. The antennas are placed half 
wavelengths apart so maximum transmission is in the forward 
and backward directions. A) Simple broadside array. 

B) Binomial broadside array. 

BASIC program to calGulate interference patterns. The 
positions XC 3 j ) of the N wave sources which emit 
waves of wavelength L, relative amplitudes A[ and phases 
PC 3 are READ in from DATA. The fundamental loop adding the 
waves from the sources is lines 230 through 290. 
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100 REM SIMPLE VERSION OF I VS ANGLE 

1 10 READ N>L 

120 F0H 1=1 TO N 

130 READ KC I ># V(I J, ACI PCI ) 

140 NEXT I 
150 DATA 2, 1 

160 DATA .5,0#1j0^ -.5# 0,1,0 
170 PRINT ••ANGLE", "INTENSI TY" 

180 FOR T=0 TO 6.28318 STEP 3.14159/18 
190 LET X0= 1O*C0S< T) 

200 LET Y0=10*SINCT) 

2 1 0 LET I 0=0 

220 FOR P0=0 TO 6 . 283 1 8 -3 . 1 4 1 5S' /8 STEP 3.14159/8 

230 LET A=0 

240 FOR 1=1 T0 N 

250 LET X=X0-XC1 ) 

260 LET Y=YO-YCI ) 

270 LET R=SQHCX*X+Y*Y) 

280 LET A=A+ACI 5^SINC6.28318*R/L + P0> 

290 NEXT I 

300 LET IO=IO+A^A 

310 NEXT PO 

320 PRINT 180*T/3. 14159,10/8 
330 NEXT T 
340 END 
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INTRODUCTION 

Past publications concerning the uses o£ the computer in field mappings 
have stressed particular geometries for the charge distributions^. This note 
presents general methods to find and map field patterns for general distributions 
of charges and current loops. Two fundamental programs are presented with il- 
lustrations of the student’s uses of both. Both prograirs use subroutines as 
building blocks so that the student can write programs using very simple calls 
to subroutines. The programs are given in appendices in the language BASIC 
as implemented on Dartmouth’s time-sharing computer system. 
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GENERAL ELECTROSTATIC FIELDS 

Th« first program allows students to map electric field lines and equi- 
potentials for an arbitrary number of point charges placed anywhere in the 
plane. The program could also map fields of line charge distributions or any 
combination of line charge and point charge distributions. The program just 
needs an expression for the electric field E (due to the charge distribution) 
everywhere in space. Thus the program can perform the calculations for uni- 
formly charged spheres, i'f one wishes 

The strategy of the program is as follows; 

The program to follow electric field lines first chooses the start- 
ing point for the field line; second, calculates the electric field in the 
X and Y directions at that point; third, uses the differential equation for a 
field line to take a small step (along the field line) to a new point in space 
Finally, the program goes back to the second step and repeats the calculation. 

For equipotontials, the program first chooses the starting coordinates 
for the equipotential; then, second, calculates the electric field at that 
point; then, third, moves a small step perpendicular to the electric field f 
to find a new point on the equipotential. Finally the calculation goes Back 
to the second step and repeats itself; 



The method only needs to calculate the electric field E at each point in 
space. The steps taken along the field lines or equipotential surfaces must 
be small enough to be treated as differentials. Ihis latter demand is more 
simply stated by saying that the numerical approximation to the integration 
of the differential equation must converge. The initial choice of starting 
coordinaire may be a set of inputted coordinates X, Y or could be the result 



of a hunting routine to find, say, some point on the 2 volt equipotential. 

Figures 1 thru 4 show the results of these calculations for several con- 
figurations of point charges. Figure 1 is the field lines and equipotentials 
three, equal, positive charges placed on the corners of an equilateral 

ERIC 



31 



\ 



- 3 - 



triangle. Notice that the starting points for the equipotentials were choosen 
in such a way that equal steps in potential occuri^ed. Figure 2 shows the field 
lines and equipotentials for a system of four, equal, positive charges placed 
on the corners of a square. Figure 3 shows the field lines and equipotential 
surfaces in the plane of a two dimensional quadrupble(two positive and two 
equal, negative charges on opposing comers of a square). Notice that the 
starting points for the field lines were not chosen quite symmetrically , Fig- 
ure 4 shows the field lines and equipotential surfaces in the plane of a one 
dimensional quadrupoleC^-wo back-to-back dipoles). 

Charge distributions even as simple as these are not usually given in 
introductory physics courses simply because of the complication of the closed 
form solutions to these problems. Using the iterative terrhniques, outlined 
above, the problem is not only tractable but even easy for the introductory 
student to underst^d and to compute by himself. 



CRITICAL PARTS OF THE ST^TEGY 

Two separate parts of these calculations deserve a little attention. 
First, how does one calculate movement along a field line? Consider taking a 
step AS from a point X,Y at which the electric field has components E^, and 
E^. Using the fact that the triangle with sides dX, dY, and dS is similar to 
the triangle having sides E , E and /E ^ + E ^ 

& a y 2c y 

AX ^ asce //e: 2+e 2) 

- X X y 

AY = ASCE //E ^+E ^3 
y X y ^ 

Finally, then, the new coordinates X,Y are just the old coordinates X,Y plus 
AX and AY, The fact that E lines are everywhere tangent to the electric 
field i leads to these simple equations for integrating along a field line. 
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This process is easily introduced to a student* The student can then use the 
procedure as a subroutine, which is general enough to be available for any 
field mapping problem. 

The second part of the calculation which deserves a little attention is 
following an equipotential . Consider stepping the distance As perpendicular 
to the electric field E at the point Use the fact that the slope of a 

line peipandicular to a line of slope, K has, in turn, the slope -l/K, Then 
for a step AS along the perpendicular to the electric field one has 

AX = -ASfE /✓B ^+E 2}' 

^ y X y 

AY = +ASCE /i€~2+E 23 
^ x' X y - 

Finally, again the new X,Y is just the old X,Y plus AX and AY, 

Using subprograms such as these, the student can write his own programs 
to perform more complicated calculations of various types. Appendix I shows 
such a program finding electric field lines and equipotentials for an arbitrary 
point charge distribution. The program is written in the Dartmouth form of 
the language BASIC. Notice that, since this form. of the language BASIC has 
global variables, the calls to subroutines are more cumbersome then they are 
in languages such as FORTRAN. This is one of the few eases in which FORTRAN 
is easier for the student to use than is the language BASIC. Some newer 
forms of BASIC will use local variables in subroutines. 



GENERAL MAGNETOSTATIC FIELDS 

The second fundamental program reported in this note calculates the 
magnetic field ^ at any point in space due to a current loop (which is par- 
allel to the X,Y pla^ne at some point Z) . TTie program simply uses the Biot- 
Savart law and integrates around the current loop. The student then uses 
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this program (a^s a subroutine) to calculate magnetic fields due to such 
current loop configuration as Helmholtz pair, or even a short solenoid (treated 
as a number of separate loops placed next to each other along the Z axis). 
Another interesting program calculates the field at some point X,Y,Z due to 
a straight wire segment. Wlien used as a subroutine, this straight wire pro 
gram can do another set of interesting current geometries. 

After calculating the magnetic field B at some point, the student can 
call the field line subroutine, considered above, and then plot field lines 
for the coil configuration. Figures 5 and 6 show plotted results from pro- 
grams such as these. Figure 5 is a magnetic field line pattern for a single 
loop around the origin. This problem is sufficiently complicated that in 
closed form it is not normally given to introductory students. Again, since 
the calculation is, first of all, iterative and, second of all, based on using 
subroutines as building blocks in programs, the introductory student has no 
difficulty understanding this approach to a fairly complicated problem. Fig- 
ure 6 is the magnetic field lines for a Helmholtz pair. Notice that the 
homogeneity in the central field region is clear, as is also the symmetry of 
the inhomogeneous part of the field. Other interesting figures can also be 
introduced. For example, if the student treats a short solenoid as a series 
of current loops placed near each other, not only does he observe (in the 
plotted field lines) the obvious homogeneity and additive character of the 
field from the various loops, but he also sees immediately the effects of 
his approximation. Because the student has made discrete what is more closely 
approximated by a continuous current loop distribution, there is imposed on 
the homogeneous field an inhomogeneous ripple with the perioda city of his 
chosen loops. 
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MAGNETIC LOOP SUBROUTINES 

The strategy of the subroutine which calculates the magnetic field at a 
point XjY,Z due to a current loop is as follows: The program, first, breaks 

the loop into segments. The program then takes the cross product of the loop 
segment and the X,Y,Z vector from the origin to the observation point. The 
program then divides by and sums the result over the loop segments. This 
subroutine is very easy to program and is also very easy for the student to 
understand. Having written the subroutine, the student simply calls this 
subroutine (and perhaps others) to plot field lines. The student uses a 
building block approach, putting blocks together to form the structure he 
wishes. Appendix II contains a loop subroutine and a program which calls 
this subroutine and calculates the magnetic field at any point in space due 
to a single loop. 



CONCLUSION 

This note has presented several programs to solve general field mapping 
problems in electrostatics and magnetostatics. Students are encouraged to 
write programs as blocks using the subroutine capabilities oif even simple 
languages such as BASIC. The student can then write driver programs (itiain 
progrMis) which are both particularly simple and also intuitively easy to 
understand. The student reponsa to systems such as this has been very good. 
Calculations based on these programs have instilled in the students a much 
deeper, and more Intuitive understanding of field patterns due to electric 
charge distributions and magnetic fields, 
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FIGURE 5 
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FIGURE 6 
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Re ferenc^ 



^J. R. Merrill American Journal o£ Physics, 837 (1969) 



Captions 
Figure 1 



Electric field lines and equipotentials in the plane of three, equal 
point charges on the comers of an equilateral triangle. The equi- 
potentials are in equal steps Ctl.2,3,4 and 5 in normalised units); 
the field lines start at equal angles around each charge. The cen- 

ter of the figure is marked with a cross. 



Figure 2 Electric field lines and equipotentials in the plane of four, equal. 
' point Charges on the corners of a square, Tlie equipotentials are in 

equal stops (tl.2.3 and 4 in normalized units); the field lines 
start symmetrically at equal angles around each charge. The center 

of the figure is the cross . 



Figure 3 



Electric field lines and equipotential surfaces in the plane of a 
t„o-dimensional quadrupole Cequal but opposite point charges on 
alternating corners of a square). The equipotentials are in equa_ 
steps except for the outermost set; the equipotentials displayed a 
i 1 5 1 0 1 5- The electric field lines do not start precise y 

symm«riUlly and their shape far from the charges is sensitive to 

the asyimnetry. 

Figure 4 Electric field lines and equipotentials in a plane containing a 

cne-dimenslonal quadrupole (two dipoles head-to-head). Again the 

-Far reeion to small asynmietries 
field lines are very sensitive in the tar r g 

in the choice of starting point. 
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Figure 5 



Figure 6 







Magnetic field lines in a central plane due to a circular current 
loop» The loop is perpendicular to the page and is centered on 
the middle; two dots mark the places where the loop cuts the page. 
The pattern is syminetric under rotations around the central field 
line * 

Magnetic field lines in a central plane of a Helmholtz pair. The 
pair is perpendicular to the page; four dots mark the points where 
the loops cut the page* The central homogeneity of the field pat- 
tern is clear as is the symmetry of the inhomogeneous part of the 
pattern. 
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EdVSUBS 



i '.>0 
1 10 
I ao 

1 30 
1 /40 

1 so 
1 60 
170 



n'KK E FI ELEv LINES & EQUIPOTENTI 

'StJl'? feldlinf.jmoa?epr:rp 

PEM c ENK > R'NY.) - C EX ^ ET > 

DF.P' -FMX C X , Y ) ==X / ( R‘4^R * R ) . 

DEF FNY CK, Y > =7'/ ( RM^R>i<R > 

PRINT "K*'j”Y“ 

REM FIELD LINE ; PART 

PRINT "INITIAL POINT <Xj.Yi: 999 



ALS USING SUBROUTINE 



TO DO KQUl-Vi ; 



I SO INPUT XliYl 

190 IF XI =999 THEN 370 



800 LET N=0 
810 LET N 1=0 
2S0 LET P‘1 = 0 



830 LET F2 = Q, 

8 AO LET FO=.l 
250 LET R=1 

S60. LET F3=FNX<X1+F1YS,Y1+Fg/S5 

P70 LET FA=FMYeXl +F1/S# Y1 +FS/a5 

880 GO SUB §i 

890 LET N1=N1+1 

300 IF NK5 THEM 340 

310 LET N1=0 ‘ 

380' PRINT XliYl 
330 LET M=N +1 
340,IFN<10 THEN S60 
350 GOTO 160 
360 STOP , 

370 REM EQUI POTENTIAL FART . 

380- PRINT "INITIAL POINT ON EQUI-V CX^YJ 999 TO END)"J 

3 90 ■INPUT XI ^Y1 

400 IF Xl=999 THEM 560 

410 REM FOLLOWS 'EOUi-V BY MOVING P E 

480 LET N=0 

4’30 LET Ni =0 

440 LET M9=. 1 

450 LET R=1 

46Q^LETlM7«FWXCXl+M5X8, Yl+M 6 /a> 

476 ; LE;T M 8 =FNYCXi +M5X2J YI+K 6/8 J 
4S0?: GOSUB #g 
:49er;^:LET . Ni=Ni,+ i 
500:- IF NT<10 THEN 540 
SlOl.LET N1=0 ' 

SSOIPRINT X1 j»Y 1 ■ , 

:S30|t.ET' M=N+ r' 

S40;^F N<10 THEN 460 
5 50i'!(i.OTO 380 
5€*.0V',Et.lD 
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'C! 






® '^S '>Fj"^^L’''^'™ESE!iT). ra=D3ST. TO MOVE TLSwL L^TiT 

- ^.. Fn,F^=n,Y nONpc; ,_nj_FT FI n rnMiy^^TTf,^ nmn^jjf 



AND 



_if |pi^' FIELD i/2 STEP AHE/^D OF PRES . FT 

-3J--XEX_£5.=5£!£.C F3* F5-i-F/a * Wa i 

56 LET Fl=Fe*F3/F5 ' * *" “~ 

-iJ — LET FF rFfl^F^/F^ ____^__ 



45 

3X 



LET 

XE12L 



oG 

sa 



RETURN 

-E ND, 



X1=X1+F1 
Yi =Y l-t-FF 
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?1QV/£H.EHP- 



S REM 
I i . R EM 



su BR.aujrjiiE 
CALCS N’EW PT 
AT CTWEM, RX 



-CXi 






....JiEXT_.X.^'L,.CLN. 
ST. M 9 ALONG 
.Yl) 



FQ, U T PCTF/j/n Al 
PERP. TO A GIVEN 



LINE 



1 S 

.zu 



REM 
„REDL 
2 S REM 

-XL =? 



3 S 

zll. 

4 S 



LET 
1 FT 



GIVE PROG.: X 1 ,Y I (PRESENT), M9=D:fST TO MOVE, AND 
fM7 = VECTOR nEETNINO GIVEN 1.J N.E AT POINT X I .Y 1- 
PROGRAM RETURNS XI, Y KNEW) AND M5 , MS = (DX , DY) 

PPnnPAtv; uses MA , Mg ,MS, M7 , ms, MQ 

M 4 = GQ R C M 7 * M 7 +M MS ) 

M.4rf'0O^f M7 /M^l ^ 



LET M5 = MS=i^ (-M8/M4') 
i FT Yt =y l+MS 



5 S LET XKXI+M 5 



SS END 



O 
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L.00 PFELD 

lOO iiEM G£dES FIELD LINES IN KZ PLANE 

1 1 o REM ARO uN D ORIGIN IN K V PLAN E--USE 

120 SUB loopb;feldlins 

130 PRINT "CUaREINT^ & KADIUS"5 

•140 'input I, a 

150 LET L0=0 

160 LET FO=.OS 

170 PKINT ’’liNITIAL POINT 

IdO INPUT Kl# 21 

190 LET N=0 

aOO LET Nl=0 

2 10 PHI NT 

aao LET KS=Kl+Fl/2 

230 LET TS^O 

240 LET 22=21+F2/2 

250 GOSUS ffl 

240 LET F3=LS 

270 LET F4=L7 

230 LET Tl=21 

290 GOSUB #2 

300 LET Z 1=T1 

310 LET N 1=N1+ 1 

32'0 IF Nl<5 THtEN 36 0 

330 LET Nl^O 

340 LET N-N+ 1 . 

3 50 PRINT K 1,7.1 

360 IF N<10 THEN 220 
370 PRINT 
330 GOTO 170 
390 END 



O 
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LOO PS 

I REM SUBROUriME — C ALC5. B OF LOOP AT Ci I v’lv'J POINT 

6 R£m give I C CO/TREN r 5 , AC liADI UGj C '-la, K'.a ) =F I eld PT-j At^n 1.0-7. OF 1.00P 

II LEM HETGHNS L S j L 6 > L 7 = C SX s BY j BK ) AT <K2 jYS^Z 25 
16 HEM LOOl^ IS PAHALLEL TO KY PLANE AT PT- E^LO 



2 1 


REM 


USES K1-4J LO-O; MO-9 




LET 


MO- lE-7=t=I 


V? 1 


LET 


K4=3. 14 159255 


36 


LET 


L1=2+K44=A/38 


^4 1 


LET 


L2=K4/ 1 5 




LET 


L5»=0 


5 1 


LET 


L6 = 0 


z>S 


LFiT 


L7 = 0 


6 1 


r'Oti 


M9=0 TO 2*K4-La STEP L2 


66 


LET 


L4=COSCMO> 


7 1 


LET 


L3 = SIKKM9 5 


7 6 


HEM 


M 1 j 2l» 3 - DL ,0N LOG P 


g 1 


LET 


M1»^-L1*L3 


S6 


LET 


M2=L l*L4 


9 1 


LET 


M3 = 0 




HEM 


L3,-9j 0 . s COORDS. OF LOOP 



tCl LET LG = Ai=L4 
•06 LET L9 = A^<L3 

1 1 1 REM Sj6 - R FROM LOOP PT . TO FIELD PT 

il6 LET MA=X8-L8 
iSl LET M5=Y2-L9' 

126 LET M6=Z2-L0 

131 LET M7 = S0PTM4:1:M4+M5*M5+M6=^HM6) 

136, LET M'3=M7*M7*M7 

141 HEM Kl^ 2, 3=CR0SSPR0DUCT 

146 LET Kl=M2*M6-M3=t'M5 ' 

ISl LET Ka=M3^M4-Ml*M6 
156 LET K3=t<5lHcM5-M2*M4 
16-1 P.EM L5,. 6^ 7=C0MPS. OF S 
166 LET L5^L5'+M0*K1/MS 
171 LET L6-L6+M0*K2/M3 
176 LET L7=tL7+M0HsK3/Mg 
181 NEXT MO 
136 HETURM 
191 END 



O 
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HSLMHOLZ 

1,00 HS:M HELMrtLOTZ PAIR B FIELD USING SUBROUTINE LOOP)! 

110 SUB LOOPS 

IfiO PRINT "CUHiiEN'1% LOOP RADI USj ■ & RIGHT SI DE COIL Z COOHD.’*; 
130 INPUT IjA, EO 

lAC PRINT "FIELD POINT OC^Y^Z)*’; 

150 INPUT xa^Y2^za 
150 LET LO=ZO 
170 QOSUB if I 
ISO LET B1sL5 
190 LET B2=L6 
aoo LET 83sL7 
210 LET LO=-ZO 
220 GOSUB if I 
230 LET B1sB 1+L5 
2 00 LET B2=B2iHL6 
250 LET B3=B3+L7 
250 PRINT . *’B FIELD:" 
aTO PKINT " C"JB15" 

280 GOTO lAO 
290 END 
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ABSTRACT 

The Schrodinger equation solutions for the sequence of potentials 
VCx) = ]x]^ for n= 1, 2, is investigated. The eigenvalues and eigen- 

functions are calculated numerically using a short program students write 
themselves. The program and its results are illustrated in the note. The 
ground state eigenvalues illustrate an interesting competition between 
kinetic energy (wavefunction curvature) and potential energy as the exponent 
n increases . 
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Introduction 

Tha Schrodinger equation for stationary states in one dimension is 
easily solved numerically using even the smallest computer. The problem 
is especially simple if the potential^ VCx)j is symmetric so that each 
eigenfunction is either of even or odd parity. The method of solution is 
very similar to iterative numerical solutions for F = ma mechanics problems. 
This note points out an interesting sequence of potentials which approximate 
an infinite square well but have rounded corners Eigenfunctions and eigen- 
values are reported for the sequence , and a simple program which students 
themselves write is presented. The sequence of lowest energies for these 
potentials demonstrate a balance between kinetic and potential energy to 
achieve the lowest total energy for the ground state. 

Potentials 

The infinite square well is a useful approximation to many physical 
problems, mostly because the problem is exactly solvable. Many physical 
situations, however, are more realistically modeled by a wall with rounded 

corners o The sequence of potentials V(x3 - |x|^ for n = 1, 2, 3, 

proaches the infinite square well but rounds the corners of the potential ^ 

Figure 1 shows several of these potentials. The curves for n = 2, 6, 
10, and 40 are shown. for x > Oj the infinite square well is shown for 
comparison. All the potentials go through the points C0»0) and (Ijl)- 
As n increases, the corner for x ^ 1 sharpens, and a larger region for 
which V(x) is essentially zero occurs. As n increases, the potential 
beyond x = 1 increases more sharply. 



Since all the potentials are symmetric around x - 0, the eigen-r 
functions have either even parity or odd parity. The ground state is even, 

and the excited states alternate odd, even, odd, . The program to intar- 

grate Schrodinger fs /equation to find an eigenfunction is very simple. The 
student demands that ^CO) * 1 (0) = ^ even states Cor - 0 and 
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(0) “ 1 fc>r odd states) a He then picks an energy and integrates the 

Schrodlnger equation - 2CVCx) - E) in a stepwise way out the x=axis. 

(All the problems discussed in this note will have^fi" = m - 1«) If E is an 
eigenvalue, the wavefunction will go to sero for large x; if not, the wave- 
function will diverge « For energies on either side of an eigenvalue, the 
wavefunction diverges to opposite signed infinities^ Hence, the student 
zeros in on an eigenvalue in three or four pairs of energies. Four figure 
accuracy of the energy eigenvalues is typical even for the simplest pro- 
grams . The wavefunctions the student derives are unnormalized. 

The programs the students write can be very shorts The basic inte- 
gration strategy is seven lines long. Figure 2A) shows a program which steps 
along the x-axis (XI) and calculates the wavefunction (PO) at each step. 

The program is in the language BASIC^ The calculation of PO uses the aver- 
age of the two first derivatives (P3 , P4j at the two ends of each interval 
Ax (LI). The calculations of the first derivatives use the average of the 
potentials (VI, V3) at the two ends of the interval. The wavefunction is 
assumed to vary slowly enough that its change over one step (of size LI) 
can be neglected* Figure 2^ shows changes in the program to achieve a more 
convergent numerical method. (One which is more accurate even for larger 
step sizes.) 

Eigenvalues - and Eigenfunctions 

Figure 3 shows the energies of the groimd state, E^, and the second 
excited state, E» (the third energy level), as functions of exponent, n, 
in the potential, V (x) = |xj . Table 1 gives numerical values for these 
energies. The sequence of energies is asympotic to the infinite square 

well C n = “) for which E = — g — in the present units. The values cal- 

m o 

culated numerically for n = 2, the harmonic oscillator, agree with the exact 

solution, E = Cm - i/2) to six figure accuracy using the higher conver- 

m 

gence program (four figures using the simple program). Figure 4 shows some 
of the wavefunctions calculated by the program. Even for an eigenvalue cor- 
rect to six figures, the eigenfunction will still diverge for large enough 
X. For eigenvalues correct to three or more figures, the wavefunction in 
the region of interest (0 ^ jx| £^1.5) Is not significantly changed by 
O „.gher eigenvalue accuracy, 
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The minimum in the values of the ground state energy, as a function 

of exponent n has an interesting physical basis e For n between 1 and about 
Sp the major effect on states of low energy (E ^ 1) includes not only the 
more steeply increasing potential beyond x - 1, but also the substantial 
broadening of the region (for x < 1) where V is essentially zero. The 
inGreasing breadth of this V 0 region allows the wavefunction to spread 
out and hence to decrease its curvature in the central area. This decrease 
in curvature is a decrease in kinetic energy, and, for small n, is larger 
than the increased potential energy experienced by the tails of the wave*- 
function in the classically forbidden region. Thus, initially as n in*- 
creases from 1, the total energy the ground state, E^, decreases. 

After n^6, most of the broadening of rhe V 0 region has been ac*- 
GOmplisb.*d and the further broadening is progressively smaller. For larger 
exponents, n, then, the cost of the increased penetration of the classically 
forbidden region rises sharply, and the total energy of the ground state 
rises accordingly* By attempting to cut off the eigenfunction tails sharply, 
the wavefunction’ s curvature and hence the kinetic energy increase again. 

For states with energies substantially greater than one, the increased 
Vv^ 0 region is of much less importance. For these states, the dominant 
effect, as h increases, is always the steepening of the potential beyond 
k = 1. For these higher energy states, the total energy should increase 
monotonically with The energy of the second excited state, E^, which 

is shown, in Figure 3, illustrates the dependence for these higher energy 
states upon the exponent n. Figure 3 also shows that the dependence of 

on n approaches that of higher energy states after the broadening effect 
of has become negligible 



Conclusion 

The sequence of potentials V (xj - 1^1^ n-1, 2, 3, --- has in- - 

teresting properties. Not only does this sequence approach the infinite 
square well (but with rounded corners), but also the stationary state en- 
ergies of this sequence illustrate tie effects of competing changes in 
kinetic energy (wave function curvature) and potential energy as the ex- 
ponent n inc^'eases. This not© has presented a simple program to find eigen- 
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functions and eigenvalues for one 
The method is easy to explain and 
in introductory quantum mechanics 
the Schrodinger equation, and has 
the physics of what is happening 



dimensional Schrodinger equation problems « 
is very general.. This use of the computer 
has removed much of the mystery of solving 
ai lowed the students to concentrate on 
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For potentials that are not symmetric, the student can be led to find 
asymptotic solutions for the wavefunctions , and then the student can in- 
tegrate in from the asymptotic region r 





TABLE 1 



Ground state energy, E^, and second excited state energy, E^, for V(x3 = |x| 
= m = 1 , 



n 


h 


^3 


1 


,809 


2.58 


2 


.707 


3.54 


3 


.675 


4,20 


4 


o 668 


4.70 


5 


.672 


5.08 


6 


■ 681 


5.40 


7 


.692 


5.66 


8 


.704 


5,89 


9 


.717 


00 

o 


10 


.729 


6,26 


12 


.753 


6.57 


14 


.776 


6.82 


16 


.796 


7.05 


18 


.814 


7,24 


20 


.831 


7.41 


30 


,898 


8.05 


40 


.942 


8.48 


CO 


1.234 


11.10 



(Ep * 6,36; B_ = 9,19) Charmonic oscillator 

with k = 2) 



CE = 16,3; E- = 29,9) 



(Eg s 30.8; = 60,5) (“ square well) 
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Figure 1. 



Figure 2 . 



Figure 3. 



Figure 4. 



Potentials V Cx) = |x|^ for n,= 2, 6 , 10, and 40. The infinite 
square well of width 2 is shown for comparison^ For n up to 
about 6, the breadth of the region where is nearly zero 

increases substantially. 

Program to solve ID, stationary state, Schrodinger equation 
problems. A) The program is in the language BASIC, The 
fundamental cal GUlations are in lines 230-^530. B) Replace- 
ment lines for a form of the program having faster convergence. 

Energy levels as a function of the exponent, n, in the potential 
is the ground state energy; E^ is the third energy level 
(the second excited state) « The minimum in near n = 6 is 
discussed in the text. Both curves are asymptotic to the in^^ 
finite square well (n - . 

Several wavefunctions for the potentials V (x) = 

Top) Ground state wavefunctions . for n = 2, 10 and 40. The 
traces are halted just as the wavefunctions start to diverge. 
Bottom) Second excited state wavefunctions for n = 2, 10 and 
40, All traces are drawn directly by the computer from data 
calculated by the program of Figure 2. 



56 



FIGURE 1 




FIGURE 



A5 SIMPLE PROCEDURE TO SOLVE ID SOHR0DIN0ER EON 

100 DEF FNVCX>=XtN9 
110 PRINT "N"l 

IgO INPUT Nf ; 

130 PRINT *'£ND OF INTERV^# STEP SIZE'*J 
140 INPUT L^LI : 

ISO PRINT ••INITIAL WAVErNt# OERIV."l 
lAO INPUT P8#P? 

ITO PRINT »'ENERay»'l 
ISO INPUT E 
IfO LET P0-P8 
eOO LET P4«P9 
SIO LET X1«0 
ego LET VlaFNVCO} 

030 LET V3-V1 

840 LET VI»FNV<XI+Li> 

870 LET Pga8«<CVl4>V3>/g - E>*^0 
8S0 LET P3-P4 

310 LET P4sP44>P8iIiL1 i 

380 LET. P0«P0+LI*€P4*P3)/8 

330 LET XI «X1 ♦LI 

340 LET NsN^Ll 

3S0I1F N«L/10 THEN 830 

3S0 LET N«0 

370 ;PR1NT XllPd 

380 ilF Xi<L THEN 830 

390 ieotO 170 

400 |£ND 



89 ^0R FASTER CONVERGENCE 

■' y- ■' - ■ ' ■ ' : 

840 LET VI*FNVCXI+L1J 

8S0 LET V4SFNVCXI+LI/4J 

260 ;LET V8-PNVCX1+LI/89 

870 LET P8s2*CC<VI^V3>/8+V8J/8 E>*PO 

280 LET P3bP4 i i 

890 LET P5«2*CC€V8+V3J/9+V49/8 -i E14P0 

300 let P64P4+PS4L 1/8 

SlO^LEf P4«P4+P2*Li 

380 LET P0iP0^LT*CCP4^P3>/ 
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This paper presents some of the preliminary results of Project COEXIST at 
Dartmouth College* Project COEXIST is an NSF sponsored project to investigate 
ways to use computers in introductory physics and mathematics teaching * The 
computer opens up areas of study not previously available at introductory levels. 
Students use the computer in a nianber of ways on homework, on individual projects, 
and in the laboratory. Students, v.hanever possible, write their own programs . On 
occasion, it is more useful to have the students programs written and saved by the 
teacher. 

Project COEXIST has made wide use of computer-connected graphic display devices 
such as X-Y plotters and cathode ray terminals. This paper presents illustrations 
of the Project'S uses of these devices. Illustrations from a broad range of study 
areas are introduced to show the versatility of these display devices. 



EQUIPMENT USED 
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X'Y plotters of the type used are made by several companies. The>- utilize 
a regular X-Y analog recorder with an interface between the recordar and the com- 
puter. The interface converts digital characters coming from the computer over 
a telephone line into voltage positions on an X-Y page. The computer software 
converts calculated values of X-Y coordinates to sets of ASCII characters. The 
ASCII characters are the elements actually sent over the telephone line. Such 
an analog X^-Y recorder is corsiderably faster in real time (and in CPU time) than, 
for example, a Cal Comp incremental system. The most important attribute of X-Y 
plotter systems and graphic systems, in general, is that the graphic or plotted 
output. is usually much more useful than tabulated numbers. While graphic displays 
do not have the resolution of tabulated numbers, they nevertheless convey a great 
deal more infomation in a short time. The cost on one of these inexpensive X-Y 
plotter systems is about $3300, of which the cost of the interface itself is about 
$2000. A typical X-Y plotter system, connected to a teletype ^d to an acoustical 
coupler is shown in figure 1. 
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FIGURE i 



Teletype and X-Y plotter as used with the time-shared Dartinouth GE-635 system. 
Software saved in the computer can be called froin any program to position the 
X-Y pen anywhere on the page. The final output can then be a hard copy paper 
plot, graph or picture. 
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A second graphic display device is the cathode ray Tektronix 4002 terniinal . 
This terminal is similar to an X-Y plotter except that it incorporates the key- 
board. Also 5 hard copy can not be obtained from the cathode ray terminal with- 
out photographic techniques or a hardware attachment. However, the cathode ray 
terminal is often very much faster than an X-Y plotter, which makes the terminal 
useful for debugging, and for those situations where the student must observe 
large numbers of curves. Since no inexpensive hard copy is available from the 
CRT display device, the CRT will not be emphasized. Nonethaless, everything in 
this paper has been performed on the CRT device. The Tektronix 4002 cathode ray 
system is shown in figure 2. 

ILLUSTRATIONS 

Illustrations from classical mechanics, fluid fields, electrostatic fields, 
geometrical optics, and physical optics have been chosen. The ways in which graphic 
display devices have been used, in as broad a range of introductory physics topics 
as possible, are shown. In classical mechanics, students have written and used 
programs on such areas as trajectory motion, strobe photograph labs (including 
non-closed integrable force laws), Keplerian and non-Keplerian orbits, relativ- 
istic. dynamics (Including relativistic motion of charged particles), acceler- 
ator simulation, and a number of scattering situationSt The following are results 
of one sophmore- level scattering simulation laboratory. 

Figure 3 shows the results of the classical scattering of point positrons off 
a model of S-state hydrogen. The atom is modeled as a point nucelus surrounded 
by a uniform negatively charged sphere. The total negative charge in the sphere 
exactly cancels the total charge of the point nucleus. The scattered particle is 
repelled from the nucleus. The figure shows the trajectories of positrons with 
an energy of ,25 of the ionization energy (13.7 volts) and for various impact par- 
ametarSo The - scattering is entirely classical . The student becomes familiar 
with the concepts of impact parameter, angular momentum, differential cross section, 
total cross section, and effective potential, in a classical system, 'fhe student 
than finds these concepts much easier to understand in quantim mechanical cases. 
Figures such as figure 3 are used in a laboratory simulation experiment. The 
student starts with hard spheres. After plotting and understanding the trajec- 
torias of hard .sphere scattering, the student than plots trajectories such as 
these shown in the figure. The student measures the angle of deviation and then 

the number of particles scattered into 20® intervals of scattering angle 
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FIGURE 2 



Cathode Ray Display terminal as used with the time^shared Dartmouth GE-635 
system. Software saved in the computer can be called from any program to 
position the light spot anywhere on the screen. The keyboard is part of the 
display terminal and letters are ’’printed’’ by the light spot on the screen. 
Hard copy Is only available by photographic techniques or from a hardware 
attachment. 
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FIGURE 3 



Classical scattering of a point positron off a model of an S=state hydrogen atom. 
The model is a point nucleus surrounded by a uniformly oppositely charged cloud. 
The whole structure is neutral. The units are normalized so that the circle is 
the out<;ide of the cloud and is one spatial unit; energy is in units of 13.7 ev. 
The student plots many such trajectories and measures angles of deviation, b is 
the impact parameter. The program illustrates a classical mechanics use of dis- 
play devices. 
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as a function of the angle of deviation* The plot is normalized in such a way 
as to make it a plot of differential cross section ^ The student gets results 
like those shown in figure 4. This figure shows the number of particles scattered 
into a given angJe as a function of scattering angle* Notice the large number of 
small angle scatterings « Since the potential is cut off by the electron screening, 
the number of small angle scatterings is large but not infinite. The student 
gets even batter data than this shown by using smaller bin sizes and more trajec- 
tories* The total cross section in the normalized units used in these plots should 

be approximately n since the normalized radius of the model atom is 1, and the 

2 

total cross section should be essentially HR . If the student adds up the results 
for the differential cross section, he t)’pically gets a total cross section be- 
tween 3.12 and 5*15. The data from this figure gives a value of 3.12. So, in 
general, students do at least as well as shown in the figure. 

Figure 5 shows this same S-state hydrogen model for electron scattering. 

That is, the scattered particles are attracted to the central nucleus. Notice 
the peculiar looping orbits. These orbits are correct. They are not due to a 
calculational error. The student must explain these effects. The answer is seen, 
most easily, in the effective potential. Figure 6 shows the effective potential 
for various impact par^etars. The physical situation is that of figure 5. Near 
an impact par^eter of .9 a normalized energy of .25), the negative charge 

just makes i,t over the bump. Since the radial velocity squared is proportional 
to the total energy minus the effective potantial, the particle slows down radi- 
ally while passing over the bump in Conservation of the angular momentum 

demands that the charge wind around the nucleus a nimiber of times. Ultimately 
the electron leaves the atom, but for angular momenta vary near ,9 it may take 
an arbitrarily long time to get away. 

In COEXIST, graphic displays have been very useful with various field pat- 
terns. The concept of vector fields has been Introduced by means of flow patterns 
in hydrodynamics. Figure 7 shows the pattern of \.i\e velocity field around a 
cylindrical object placed in a uniform stream. Students plot such flow patterns 
for a nimber of objects. Tlien they find densities and directions of lines to 
derive the relative velocities at different points in the pattern. In this way 
the student not only acquires some information about fluid flow, but also about 
various ways to represent vector fields. Vortex fields in fluid flow can often 
motivate curl discussions very well. 

Figure 8 shows the jilot of the electric field lines and equip otentials 
around a two dimensional quadrupole. This program allows placing point rharges 
Q inywhare in the plane. The program then follows the field lines by integrating 
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FIGURE 4 
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Differential cross section as a function of scattering angle for trajectories 
such as those of figure 3. The scattering angles are grouped in 20 degree inter- 
vals. The differential cross section is then the number in a bin times the impact 
parameter b, divided by sin0A6. The small angle scattering is large but not in- 
finite since the potential is cut-off by the charged cloud. 
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FIGURE 5 



T.S 



Classical scattering of a point electron off the model of S-state hydrogen. The 
looping orbits are real. Near an impact parameter of .9 for energy of ,25 Cnorm= 
alized units) the particle takes an arbitrarily long time to escape from the atom. 
The program illustrates another of many uses of computer graphics in classical 
mechanics . 
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Effective potential for several impact paiameters for electron scattering off 
hydrogen. The curves correspond to the trajectories shown in Fig* 5* Since the 
(radial velocity)'^ is proportional to , the particle can take a long time 

moving in or out over the bump in During this tiine^ conservatio i of angular 

momentum forces the particle to wind many times around the nucleus. 
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Velocity". l|;nes- flow around a cy-lindrycal object. Tfie student becomes used to field 
l^e concejpts by\ leastirm^ tfe CrelativeX velociti'es at various points in tbe velocity field. The 
progr.am is one of several dealing witK fluid flow and velocity fields. 
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Electric field lines and equipotential surfaces in tlie plane of' a two-dimensional quadrupole. The program allows 
placement of charges anywhere in the plane. Field lines are followed by integrating dx/E = dy/E equipoten- 
tial surfaces are followed by hunting. The program illustrates one of many ways computer^graphic^ displays are 
used to enhance the student's intuition about electrost ’;tic field problems. 




follov/s the equipotential contour. Programs have also been developed which follow 
the equipotential by moving everywhere perpendicularly to the local electric field 
lines* These programs and plots are useful in that the student very quickly gains 
an intuitive grasp of the meaning of the abstract concepts of field lines and equi- 
potential surfaces « Figures 9^ 10 and 11 show other charge distributions that 
students have found interesting and useful* 

In geometrical optics we have developed programs to illustrate tracing 
principal rays through a thin lens or spherical mirror optical system. By plotting 
rays for several systems the student understands imaging much more quickly. We 
have also developed true ray tracing programs in which the student can place any 
number of spherical interfaces between media anywhere in the plane. The student 
then starts the ray at some angle and at some position , and the program traces 
the ray through the system. This program demonstrates nicely the various forms 
of abberation, as shown in figure 12, The shift in focal point dae to spherical 
aberration is apparent * 

Figure 13 shows another application of the computer to introductory geometri- 
cal optics “ mirages. The program allows the index of refraction of the mediuin to 
be a function of height. The figure shows the results fur a model index of re- 
fraction near heated ground. The sheet of less dense air near the ground pro- 
duces a second, inverted, virtual image of any object above the surface. One 
observes two objects - one at the true position of a tree; the second a mirage 
image. This program can also be used to demonstrate *^ooming^^. Looming is an 
upside down mirage and probably explains most sightings of the "Flying Dutchman". 
Looming is due to a layer of warm air sandisiched between layers of cold air. The 
program also illustrates reflections of radio waves off the ionospheric F layer. 



Principle. The student, for example, actually shows that an N-slit diffraction 
pattern is produced by adding up circles centered on the slits. One can even 
produce a fairly good single slit pattern with relatively short times. 

Such programs have also led to more applied physical optics. Students have 
dealt with arrays of radio antennas placed to maximize the directionality for 
transmission or reception. They have studied broadside arrays and int arfarometric 
radio telescope arrays • In these programs the student compares plots of intensity 
versus angle, and intensity versus position on a screen. Figure 14 is a plot of 
intensity versus angle for three line sources placed on a equilateral triangle of 
side length = 1/2 wavelength. 



In physical optics, we have programmed a number of applications of Huygen’s 
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FIGURE 9 

Field lines and equipotential surfaces in a plane containing four equal positive point charges on the corners 
of a square. Notice that close to the charges and far frora the charges the equipotentials are circles. At 
intermediate distances the equipotentials in the plane take the general shape of the charge distribution. 
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Field lines and equipotential surfaces for two back-to-back linear dipoles. 
The complete set of surfaces are generated by rotating the figure about an 
O axis through the charges . 
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PIGUliE 12 

Rays traced tfiiougH a tfii'cE. lens. T&e effects' of aBerration axe apparent in tRe focal point, TTie program 
allows any'.nuiuEer of spherical interfaces Between media to Be placed anywfiere in tRe plane. Rays are then 
started at any position and angle to the axis and are traced through the system,. 





FIGURE 13 

Mirage effects wEen tKe index of refraction varies with height above the ground. Due to 
total internal reflection in the heated air layer, an observer sees two objects - one above 
the ground; one B'elow, The program is one of a series of various ray-tracing illustrations. 




FIGURE 14 



Intensity o£ emitted radiation versus angle from three synchronous sources 
pjaced at the corners of an equilateral triangle of side X/2* The program 
allows the student to place any number of sources anywhere in the plane 
and to choose relative intensities and phases. 
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A laboratory application of the computer in introductory physical optics 
is shown in figure IS. This figure shows Cnormalized) N-slit diffraction patterns 
from 1, 2j S and 10 slits using the geometry and wavelength given for an intro- 
ductory lab. The student measures the intensity versus position across the ob- 
servation screen of an N-slit experiment using a photoresistor and a He-Ne laser. 
He plots his results on top of a theoretical plot such as figure 15. The agree- 
ment is good. 

CONCLUSION 

This paper has illustrated a ntunber of uses for graphic display devices 
in introductory physics teaching. All the illustrations u.sed inexpensive X-Y 
plotter systems connected to the computer. The programs have also been performed 
on a CRT terminal. Some of the plots also have been performed on the teletype 
Itself. Hopefully, a large number of exMiples of the use of graphic display de- 
vices illustrates the importance of these devices better than volumes of words 
and perhaps these examples will trigger other examples in the reader’s mind. A 
glance at a plot often leads to more understanding than a long look at a series 
of numbers in tabular form. 

A short comparison of the three methods of plotting we have used might be 
helpful. The CRT device is faster than either the X-Y system or the teletype. 

It is useful for debugging and for situations where large numbers of curves must 
be displayed quickly. The CRT is. Incidentally, essentially silent when compared 
to a teletype, but also relatively expensive when compared to a teletype and an 
X-Y plotter. The X-Y recorder plotting system is much faster, in general, than a 
teletype and has much higher resolution. It automatically gives the student a 
hard copy of his results for further deduction.. It could well become the work 
horse of introductory teaching applications. The plots are easy to handle, scale 
and even label on the computer if the student wishes. The teletype itself can 
be used to plot some curves; often however, the resolution available is insuf- 
ficient to show anything but coarse behavior. Nonetheless teletype plotting is 
a useful way to introduce students to plotting in general and can be very useful 
in debugging plotting routines fairly quickly. 




D ^ ^ 

Q H 

t/j bo 




FIGURE 15 



